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We identify a generic new form of collective flavor oscillations in dense neutrino gases that amounts 
to a self-induced parametric resonance. It occurs in a homogeneous and isotropic ensemble when a 
range of neutrino modes is prepared in a different flavor than the neighboring modes with lower and 
higher energies. The flavor content of the intermediate spectral part librates relative to the other 
parts with a frequency corresponding to a typical Am^/2_E. This libration persists in the limit of an 
arbitrarily large neutrino density where one would have expected synchronized flavor oscillations. 
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PACS numbers: 14.60.Pq, 97.60.Bw 

I. INTRODUCTION 

Flavor oscillations in dense neutrino gases exhibit col- 
lective phenomena caused by neutrino-neutrino interac- 
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simplest example of a dense neutrino gas that is homo- 
geneous and isotropic, studies of collective oscillations in 
the two-flavor context amount to solving the nonlinear 
equations of motion (EOMs) for a set of flavor polariza- 
tion vectors P, , 



where 



P, = H, X P, 



H,; = + ^P . 



(1) 



(2) 



The "effective magnetic field" B is a unit vector in flavor 
space and uji = l\m? /2Ei the vacuum oscillation fre- 
quency for a neutrino mode with energy Ei. The total 
polarization vector of the ensemble is 



(3) 



Finally, a* ~ V2Gf Ui, is a typical neutrino-neutrino in- 
teraction energy. Its exact definition depends on how we 
normalize the polarization vectors, but for our purpose it 
is simply an adjustable parameter of dimension energy. 

Two main questions are of interest: What are the col- 
lective forms of motion in a dense neutrino gas? What 
happens when /i slowly decreases, mimicking the expand- 
ing universe or the neutrino density variation with dis- 
tance from an astrophysical source? We concentrate on 
the former question and identify a new mode of collective 
oscillations that so far has gone unnoticed. 

The simplest form of collective motion consists of "syn- 
chronized oscillations," meaning that modes with differ- 
ent frequencies LOi oscillate with a common frequency 
Wsync if M is sufficiently large [4, 6, 8]. The main physical 
idea is that the Hamiltonian vectors H; in Eq. ^ are 
dominated by /iP, forcing all Pi to follow P and thus to 



a common precession around B. It is this picture that 
we will see is not always complete. 

Another case is the "pure precession mode" where 
each Pj is collinear (parallel or antiparallel) to 
[13, [H, [H, [131. Even though H, depends on all P, 
through P, such a self-consistent solution exists for any 
strength of ji. All lie in the plane spanned by B 
and P. This plane rotates around B, all P; being static 
within it. For — > oo this solution requires all P^ to be 
collinear and then is identical with a synchronized solu- 
tion. Conversely, beginning with collinear P^ and slowly 
reducing /i takes us adiabatically through different pure 
precession modes. 

Finally, the "pendulum in flavor space" represents a 
class of solutions relevant for /i not too large |2j, [23| ■ 
The simplest case involves two vectors Pi, 2 that initially 
point in opposite directions. The dynamics of this system 
is equivalent to a gyroscopic pendulum. Once more, for 
/i — > 00 the motion is a common synchronized precession 
without nutation unless initially P = Pi -I- P2 = 0. 

It was always assumed that, unless initially P = 0, the 
large-/x behavior is a synchronized precession. Likewise, 
it was assumed that the pure precession mode is stable. 
However, while a small perturbation may lead to small 
oscillations around the ideal solution, it may also lead 
to an exponential deviation. Numerical studies always 
found stable behavior for the pure precession mode and 
synchronized behavior in the large-/i limit. However, we 
will see that these findings depend on special choices of 
initial conditions and are not generic, although probably 
most relevant in the supernova context. 

It is easy to see that the pure precession mode and 
its large-^ limit need not be stable. A pure precession 
mode with each P^ collinear with is a self-consistent 
exact solution. Realistic initial conditions, however, be- 
gin with all Pi collinear with each other, representing 
the assumption of all neutrinos being prepared in weak- 
interaction eigenstates. On the other hand, the initial /i 
may be large but must be finite. This initial condition 
is not an exact pure precession mode because each P.^ 
sports a small angle relative to its and thus moves 
on a precession cone with a small but nonzero opening 
angle. Therefore the total P itself can not be static in 
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the co-rotating plane. Each P; precesses around with 
an approximate frequency /i, so P itself must "vibrate" 
with a similar frequency. 

We thus have a typical situation for a parametric res- 
onance: Each Pi precesses around with the approx- 
imate frequency fi and Hi itself vibrates with a similar 
frequency. Therefore, while the precession cone of each 
Pi follows its Hi, it is not assured that the opening an- 
gle stays small. The exact outcome depends on how the 
evolution of the P^ feeds back on P and thus on H^. 

To identify the simplest system showing such behavior 
we return to the EOMs of Eq. ^ and note that they 
have two exact invariants. One is the angular momentum 
along B, so B • P is conserved The other is the 

energy [H, [2^ 

B • M + ^ P^ (4) 
where the "total magnetic moment" is 

M = ^c^,P,. (5) 

i 

When ^ —^ 00 energy conservation implies that P^ and 
thus I P I are conserved [22| . P then precesses around B as 
a collective object, representing the usual synchronized 
oscillations. 

This does not imply, however, that all P^ remain 
coUinear with P. There can be internal motions among 
them such that |P| is conserved and P still precesses as 
a collective object, yet the individual P^ move relative to 
each other and relative to P. 

For two polarization vectors no internal motion is pos- 
sible without modifying |P|. Moreover, since |Pi| is con- 
served, the motion of each P^ is described by two angles, 
so we have a total of four degrees of freedom. With two 
exact constants of the motion there remain two degrees 
of freedom, representing, for example, the "nutation an- 
gle" and the "precession angle" of the flavor pendulum. 
There is no room for new forms of motion. 

If we have three or more polarization vectors that are 
initially aligned, nothing new happens either because the 
approximate conservation of |P| in the large-// limit re- 
quires that the Pi remain aligned. 

Thus we need three polarization vectors (oscillation 
frequencies wi < ^2 < 1^3), one of them initially anti- 
aligned with the others. If the "flipped" polarization vec- 
tor is number 1 or 3 we are back to the flavor pendulum 
because in the large-/x limit two neighboring aligned vec- 
tors act roughly as one average mode. This is not possible 
if the flipped vector is number 2, so the relevant initial 
configuration is either up-down-up or down-up-down. 

The main purpose of this paper is to study the most ba- 
sic system of this sort and some simple generalizations. 
In Sec. In] we solve analytically a very symmetric sys- 
tem consisting of three modes. In Sec. IIIII we consider 
its counterpart for a continuous spectrum of modes and 
solve it analytically. We discuss the implications of our 
findings in Sec. IIVI 



II. THREE POLARIZATION VECTORS 

A. Basic configuration 

The basic configuration showing a self-induced para- 
metric resonance consists of three polarization vectors 
with frequencies uji < UJ2 < For simplicity we assume 
equal frequency spacings j = uj^ — uj2 — UJ2 — ^i- The 
common precession of the system around B is irrelevant, 
so we can choose UJ2 = 0, essentially going to a rotating 
frame around B. The polarization vectors are P±,o with 
oscillation frequencies ±7 and and the EOMs are 

P± = (±7B + fiP) X P± , 

Po - mP X Po • (6) 

Moreover, we assume that all polarization vectors have 
unit length and that initially P± point in the positive, 
Po in the negative z-direction. 



B. Numerical examples 

We first illustrate the behavior of this system with a 
numerical example. We take B to be tilted relative to 
the z-axis by an angle 29 with cos 26* = 0.5, 9 itself play- 
ing the role of the neutrino mixing angle. In Fig. [T] we 
show the evolution of the projections of P± and Pq on 
P for the interaction strength fj, — 10 7. The vector Pq 
evolves from its initial anti-alignment with P to complete 
alignment and back, and so forth periodically. The other 
two vectors evolve similar to each other such that B • P is 
strictly conserved and the length |P| (uppermost curve) 
is approximately conserved. 

The interesting internal evolution is relative to P, so a 
large complicates the system by projection effects. We 
may use the opposite extreme of a vanishing 6 because 
the existence of the new collective mode can not depend 
on this choice. A vanishing 6 has the only disadvantage 
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Time [271/7] 

FIG. 1; Evolution of three polarization vectors for jj, — 10 7. 
We show P± and Po, projected on the direction of the total 
P. The uppermost curve is jFj. 
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FIG. 2: Evolution of three polarization vectors as in Fig. [T] 
now with a small mixing angle and = 3 7 (upper panel) and 
fi — 10 "f (lower panel). The dotted line corresponds to the 
interaction energy (see text). 



that the system can not start moving if the polarization 
vectors begin exactly collinear with B. However, the 
motion can be excited by using a small but nonzero 9 
or by some other disturbance. 

In Fig. [2] we show the evolution for an extremely small 
but nonzero mixing angle using /i = 87 (upper panel) 
and fj, = IO7 (lower panel). We do not show the long 
quasi-static phase where the system starts moving. One 
difference to the previous case is the absence of any "wig- 
gles" even for a moderate /u. For larger /j, the variation 
of the interaction energy (/x/2)P^ becomes smaller. We 
show as a dotted line the quantity (/i/2) (P^ — 1). For 
II ^ 00 the two energy components B • M and (/i/2) 
seem to be separately conserved. It is striking that the 
"flipping time scale" remains the same. 

We sketch the motion of the polarization vectors in 
Fig. [31 The vector Pq oscillates relative to P, its relative 
orientation being described by the angle tp relative to the 
negative P direction. If we use coordinates where the 
x-direction is defined by Pq we can describe the motion 
by the angle ip alone in the form 



= 



and 




(7) 



where s — simp and c = cos (p. In this way we always 
have P = (0,0, 1). The energy components B • M and 
(/i/2)P^ do not depend on (p. 




FIG. 3: Internal motion of three polarization vectors. 



C. Analytic solution 

Equation ([7]) must be a good approximation to the so- 
lutions of the EOMs, so all that is missing is the time 
evolution of (p. To find it we simplify the EOMs further 
and observe that the component of P along B is strictly 
conserved. Our choice of unit length for all vectors there- 
fore implies P = B -I- Pj_ where P_l is the P component 
transverse to B. Therefore, the EOMs are 



P± = mP±xP± + (m±7)BxP±, 
Po = AiP± X Po + mB X Po . 



(8) 



We observe that there is a common precession with fre- 
quency ^ around B that represents the fast internal pre- 
cession of all Pi around P. We can take out this common 
motion by going once more to a rotating frame so that 
finally the EOMs for the internal motion are 



P± = nP± X P± ± 7B X P± , 
Po = mP± X Po . 



(9) 



We have achieved that fast precessions that are com- 
mon to all polarization vectors have disappeared. The 
polarization vectors see the usual B, leading to normal 
precessions, and in addition a transverse effective field 
represented by Pj_ that tilts them. 

This is analogous to the textbook example of mag- 
netic resonance where an electron spin precesses fast in a 
strong external B field. Applying even the tiniest trans- 
verse B field that rotates with the original precession 
frequency leads to a situation where in the co-rotating 
frame the spin only feels the transverse component and 
completely reverses relative to the direction of the large 
external field. 

In this example the rotation speed of the transverse 
B component must be adjusted to match the electron's 
precession frequency. In our case this resonance condi- 
tion arises automatically by the nonlinear nature of the 
problem where the transverse field is provided by the 
precessing polarization vectors themselves. 

To solve the EOMs it is useful to introduce the differ- 
ence vector 

D = P+-P_ (10) 
and study the motion of the three vectors Pq, D and P^. 
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Their EOMs are 

Po = /iP± X Po, (11) 

b = mP± X D + 7B X (Px -Po), (12) 

Pi = 7BX fo-^Pi) . (13) 



The vectors Pq and D have lengths of order unity, 
whereas in the large-/z Hmit P_l must be small. Energy 
conservation implies that it must be of order (7//^)^^^ 
or smaller. However, since P^ appears on both the l.h.s. 
and r.h.s. of Eq. (|13p and since it appears multiplied with 
fi on the r.h.s., both sides of this equation are of the same 
order in 7/yU only if to leading order 

Pi = -D. (14) 

Here we have assumed that to leading order D has no 
component along B because P-|_ and P_ are expected 
to behave symmetrically. Inserting this into Eqs. (jlip 
and (|12p and neglecting terms of order 7//! we find 

Po = 7DXP0, 

b = -7B X Po. (15) 

To leading order in 7//^ we have found a closed system of 
EOMs for Po and D alone. The strength of the neutrino- 
neutrino interaction has disappeared, only 7 appears as 
an energy scale. 

In our coordinate system Po moves in the x-z plane 
whereas D is along the y-direction. Describing the ori- 
entation of Po with the angle (p as in Eq. ([7]) , the EOMs 
finally become 

ip = ■jD, 

D ~ —7 sin 93. (16) 

We recognize that (p and D play the role of canonically 
conjugate variables and that these EOMs follow from a 
Hamiltonian 

Hiip,D)=j[^D^ + {cosip^l)] (17) 

by virtue (p — dH/dD and D — —dH/dp. Therefore, 
we have the EOM of a pendulum with frequency 7. 

The initial condition P± aligned with B and Pq anti- 
aligned corresponds to the initial condition D = and 
if = 0. In other words, the pendulum begins in an in- 
verted position and obeys 

Y=7'(l-cos^). (18) 

This agrees well with numerical examples. 

Using the vectors of Eq. ([7]) with any value ipo as an 
initial condition amounts to Dq = [2(1 — cos^^o)]^^^ or 
(^0=7 [2(1 — cos(fQ)Y^^. The system follows the same 
solution (up to corrections of order 7/^) for any <po, ex- 
cept that the pendular motion begins in some phase of 



the oscillation other than the inverted position. The pen- 
dulum begins with an initial excursion ipo and an initial 
velocity (po, both adjusted such that it will reach the in- 
verted position. 

We can also choose other initial conditions with an 
arbitrary ipo and no initial velocity {Dq = 0) . In this case 
we have a pendulum that oscillates between two maximal 
excursion angles without reaching the inverted position. 
In terms of the polarization vectors this corresponds to 
P-t being initially aligned with B whereas Pq begins with 
a nonvanishing angle ipQ. 

We can choose this initial angle such that all three po- 
larization vectors are initially almost aligned. In this case 
we obtain small-amplitude oscillations. This shows that 
our solution is also relevant for three initially aligned po- 
larization vectors, the case where one finds synchronized 
oscillations. In the limit /x — > 00 and exactly aligned 
vectors nothing happens. However, for a finite fi and a 
small disturbance these vectors do oscillate. Here the 
small opening angles of the precession cones caused by 
the imperfect initial alignment between P^ and also 
cause a parametric resonance, but in the direction of the 
small opening angles becoming smaller, then back to the 
original angle, and so forth. So the parametric resonance 
always occurs, but has a visible effect only if Po is ini- 
tially anti-aligned with P±. 

We finally remark that the initial orientation of the 
three polarization vectors relative to B is irrelevant. 
The initial conditions up-down-up or down-up-down are 
equivalent. In contrast to the gyroscopic flavor pendulum 
the neutrino mass hierarchy is here irrelevant. 



D. Large mixing angle 

We return to the case where P± and Po are initially 
at some large angle 29 relative to B, corresponding to a 
large neutrino mixing angle 9. Armed with the insights 
gained in the previous section we note that we should 
study this situation in a coordinate system co-moving 
with P and co- rotating with Po around P. Now the 
vector B rotates with the large frequency fi around P, its 
transverse component averaging to zero. In other words, 
we may consider the rotation-averaged EOMs in the same 
spirit as for the ordinary matter effect [2^, ^4^ . 

The discussion of the previous section remains un- 
changed except that only the component of B along P 
contributes. This amounts to using the effective oscil- 
lation frequency 7 cos 29. One can verify in numerical 
examples that this is indeed what happens. Even com- 
paring our Figs. [Hand [2] reveals to the naked eye that a 
large 9 value increases the oscillation period. 

We conclude that the essential dynamics of our system 
does not depend on 9. For our theoretical study it is most 
convenient to use 9 = whereas in numerical examples 
one uses a nonvanishing value that provides the necessary 
initial disturbance to start the motion. 
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E. Matter 

Ordinary matter has the same effect on aU modes 
and therefore can be removed by going into a rotating 
frame (2^ . Since we anyway only study the internal mo- 
tion of the Pi and since we already go to a rotating frame 
to achieve this simplification, nothing new happens by 
the extra matter-induced rotation. Therefore, the effects 
discussed here are not modified by the presence of dense 
matter except for details of the initial disturbance caused 
by the fast-rotating B. 

F. Antineutrinos 

Antineutrinos are most easily included in the EOMs as 
modes with negative frequencies cu = — |Am^/2£'| 28], so 
in the most general case both neutrinos and antineutrinos 
are present. On the other hand, we constantly switch 
between rotating coordinate systems, shifting the zero- 
point of frequency at convenience. After each shift the 
behavior of the abstract system is the same, whereas its 
physical interpretation changes. 

The polarization vectors used here are similar to the 
neutrino flavor isospin vectors of Duan et al. [221] in that 
for an antineutrino (negative frequency) "spin up" means 

if for a neutrino (positive frequency) "spin up" means 
i>e- Shifting between rotating coordinate systems there- 
fore changes, for example, a Dfj. to a i/g. 

We repeat, however, that the interpretation of the po- 
larization vectors is irrelevant with regard to the dynam- 
ics of the system. It is easiest to think of all modes as 
representing neutrinos in that any given distribution can 
be thought of as stemming from a neutrino distribution, 
suitably shifted to a rotating frame. Therefore, in an 
abstract discussion of the EOMs and their solutions, an- 
tineutrinos need not appear explicitly. 

G. Fewer symmetries 

We have studied a very symmetric system, but the gen- 
eral behavior persists if we vary the relative lengths of the 
vectors P±,o and the frequency splittings. However, for 
any motion of this sort to be possible the three vectors 
must be able to move relative to each other while con- 
serving P. If P is conserved, energy conservation implies 
that we also need to conserve B • M. 

Once the lengths of the three vectors are specified, 
their motion is described by two angles each, so a total of 
six angles. Our conditions provide four constraints, leav- 
ing two degrees of freedom: an irrelevant overall preces- 
sion angle around P and one angle describing the internal 
configuration. 

The general conditions for a solution are not particu- 
larly illuminating, but we remark that, if P+ and P_ 
have equal lengths, the parametric resonance requires 
|Po| < 2 |P±|. This is verified in numerical examples. 



III. SPECTRUM OF FREQUENCIES 
A. Numerical Example 

Next we consider a broad spectrum of modes. For con- 
venience we scale all frequencies with a scale 7 and use 
as a time coordinate, i.e., we use dimensionless frequen- 
cies. Without loss of generality we consider the range 
— 1 < uj < +1. We study the continuous generaliza- 
tion of the previous three-vector example, using an ini- 
tial spectrum that is of box shape and has a flipped part 
in the middle (shaded spectrum in Fig. |4]). We let this 
ensemble evolve with a fixed large /x, using a B that is 
collinear with the polarization vectors except for a small 
mixing angle that triggers the motion. All modes oscil- 
late with the same phase and reach their maximum excur- 
sion at the same time. We show the maximum excursion 
spectrum as a thick line in Fig. 31 The result is clearer 
when we "fiip" the central part of the spectrum (dashed 
line). Within numerical accuracy, this compound curve 
is a Lorentzian resonance curve. 

Therefore, the different parts of the neutrino spectrum 
librate relative to each other. In particular, the neu- 
trinos with intermediate frequencies oscillate completely 
between two flavors with a typical vacuum oscillation fre- 
quency. The neutrino density does not enter as long as 
it is large in the sense /i S> 7. 

We can make the "flipped" part of the spectrum nar- 
rower, obtaining a narrower Lorentzian. However, we 
can not make it arbitrarily broad in the same way as in 
the three-vector example the middle vector could not be 
longer than twice the length of the peripheral vectors. 
Here we find that the width of the flipped part must not 
exceed half the total box width. The transition is sharp in 
the sense that if the flipped part is slightly broader than 
half the total, nothing happens, if it is slightly smaller 
we find complete reversals and back of the central part 
of the spectrum. 
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FIG. 4: Initial spectrum (thin line and shaded) and spectrum 
of maximum excursion (thick solid line), assuming a large 
fixed fi. The dashed line is the "flipped" version of the thick 
spectrum. 
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B. Analytic solution 

To determine the exact solution in the large-/i Umit we 
note that the EOMs in analogy to Eq. ^ are 



/iP_L X + X 



(19) 



where each polarization vector is characterized by its fre- 
quency Lu. We have assumed again that the mixing angle 
is small and therefore B collinear with P — J dui P^, ■ 
Integrating both sides over / div leads to 



P_L = B X (M-^PP_l) 



(20) 



where P — \P\ is here not assumed to equal unity. The 
"magnetic moment" is M = J duj lo Pui ■ The same ar- 
gument as in the three-vector example implies P_l — 
M/(/iP) to lowest order in We use spectra that 

are symmetric relative to a; = so that M to lowest 
order has no component along B. So the EOMs are 



cjB 



M 

"P 



X P„ 



(21) 



Once more we have eliminated the large frequency /x. 

In this form it is clear that we expect a resonance shape 
of the oscillation pattern as a function of oj. Only the 
mode with w = is exactly on resonance in that the 
transverse B field (here M) is exactly co- rotating with it, 
or in our co-rotating coordinate system, does not move 
at all. The other modes precess around B relative to the 
M direction. They are not exactly on resonance. 

For an explicit example we return to a box spectrum 
like Fig. |4] and define the function 



+ 1 for /3 < \uj\ < 1, 
-1 for < |cj| < (3. 



The initial polarization vectors are 



(22) 



(23) 



and so P = |P| = /dcj = 2 (1-2/3). We further define 
the Lorentzian function 




1 



(c^/r)2 + 1 

and fix its width by the condition 



which implies 



r = 



(24) 



(25) 



(26) 



This is only possible for /3 < i, precisely the case where 
we find the parametric resonance. 



As in the three-vector example we use a coordinate 
system where M is oriented along the y direction and 
the oscillation of Po takes place in the x-z plane. Once 
more we describe this vector by an angle such that 



sm ip 


- cos (/? , 



(27) 



Inspired by the numerical example and by the symmetries 
of our system we guess 




sm ip 



(c^/r)y2(r~cos</,) 
1 — cos (f 



(28) 



Note that Eq. (^5]) and the symmetry of /^j and s^j guar- 
antee that P = J duj P^ is the same for any tp. 

Next we integrate Eq. (f2T|) over J duj uj so that on the 
r.h.s. the term M x M drops out and we find 



M = B X dLOLu-^P. 



(29) 



The vector M has only a y-component that is found to 
be rP-\/2(l — cos (p). The vector on the r.h.s. also has 
only a j/-component F^Psint/j. Therefore, once more the 
excursion angle evolves like an inverted pendulum. 



y/2{l- cos ip). 



(30) 



28l) indeed solves the 



It is now easy to verify that Eq. 
EOMs of Eq.' 



C. Single-crossed spectrum 

The behavior found in the previous section is similar 
for less symmetric arrangements. However, the "flipped" 
part of the spectrum must not be too broad and must 
not get too close to the edges of the overall box. In 
other words, the two "wings" of the spectrum must be 
large enough to support the parametric resonance, but 
we have not worked out the general condition. 

The crucial feature is to have a "double-crossed" spec- 
trum s{lo), representing the initial z-component of P^;. 
The spectrum must first cross from positive to negative 
and then back to positive values or the other way round. 
One can also construct multiple-crossed spectra and finds 
more complicated parametric resonance patterns. 

On the other hand, this form of collective motion can 
not occur for a single-crossed spectrum where s{lo) < 
for LU < ojcross and s{uj) > for lo > Wcioss- Here B • M 
is maximal, so moving any polarization vector makes 
it smaller. This is easily seen if we go to a rotating 
frame where Wcross = 0. This has the effect B • M = 

B • WiPi ^ B • (t^i - Wcross)Pi = B • M - WcrossB • P 

and since B • P is a constant of the motion, we have only 
added a constant term to the energy. In the new system 
with Wcioss = all modes with negative frequencies get 
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multiplied with a negative s{lo), the positive ones with a 
positive s{uj). Therefore, any motion of P^^ lowers B • M. 
The opposite is true if the spectrum crosses from posi- 
tive to negative values or if we take B to point in the 
negative z-direction. In other words, even though for a 
single-crossed spectrum one could move the polarization 
vectors such that |P| and (/i/2)P^ remain unchanged, 
the quantity B • M is extremal. 

A single-crossed spectrum is the prototype for a gyro- 
scopic flavor pendulum that is driven by the exchange of 
energy between (/i/2) P'^ (playing the role of kinetic en- 
ergy) and B • M (potential energy). In the large-^ limit 
the oscillations are synchronized and, if the mixing angle 
is small, nothing visible happens. 

If one assumes that a source only emits neutrinos and 
antineutrinos of one flavor, e.g. a flux of Ve (positive lo) 
together with a flux of Vf, (negative uj) one has a single- 
crossed spectrum in that the polarization vectors for j/g 
are "spin up," those for "spin down" as explained in 
Sec. IIIFI Therefore, in this situation the parametric res- 
onance does not play any role. 



D. Nonisotropic system 

The parametric resonance is probably not possible if 
the neutrino ensemble is not isotropic. Assuming axial 
symmetry around some direction we can describe every 
mode by its frequency u and the velocity projection v on 
the symmetry axis. The EOMs are in this case [1^ [25j 

where 

H„,„ =ujB + ^l{P - vF) (32) 

and 

P = JdujdvP^^y, F = JdujdvvP^^y. (33) 

Even if the vectors P and F are collinear, there is no 
longer a co-rotating frame where all polarization vectors 
would see essentially the same effective magnetic field. 
If one of them were on resonance with regard to some 
transverse field, the others would precess with different 
angular velocities where the differences are much larger 
than order uj in the large-// limit. 

We have studied a few numerical examples where we 
assumed a homogeneous ensemble with a "half-isotropic" 
distribution, i.e., all modes with velocity components in 
the positive direction of the symmetry axis were isotrop- 
ically occupied, those with negative velocity components 
were taken to be empty. The parametric resonance in- 
deed disappeared for all tested examples where it oc- 
curred in the corresponding isotropic case. 



IV. DISCUSSION 

We have shown that in a dense neutrino gas, neutrino- 
neutrino interactions can lead to a parametric resonance 
that causes different parts of the flavor spectrum to os- 
cillate relative to each other. This effect persists in the 
large-density limit where one would have expected syn- 
chronized oscillations. The new pattern is an internal mo- 
tion among the ensemble of polarization vectors that still 
precess together as a collective object, but do not retain a 
common orientation relative to each other. The strength 
of the neutrino-neutrino interaction docs not appear in 
the oscillation frequency that is determined entirely by 
the spectrum of vacuum oscillation frequencies. 

However, the parametric resonance has a macroscopic 
impact only if the initial flavor spectrum s{lo) is at least 
"double crossed" in that there must be a spectral fla- 
vor sequence of the form up-down-up or the other way 
round. We use polarization vectors similar to "neutrino 
flavor isospin vectors" where, for example, a 9^ mode is 
"spin down" if v^. is defined as "spin up." Therefore, 
a neutrino flux initially consisting of neutrinos and an- 
tineutrinos of a single flavor represents a single-crossed 
case: s{uj) changes sign once at a; = going from antineu- 
trinos to neutrinos. Most numerical studies of supernova 
neutrino oscillations used such single-crossed spectra and 
the parametric resonance did not show up. 

A dense gas of neutrinos and antineutrinos of different 
flavors in kinetic (but not chemical) equilibrium is also 
single crossed. For a Fermi-Dirac distribution the spec- 
tral shape and amplitude are fixed by the temperature 
and chemical potential. If the distributions for two fla- 
vors are described by the same T but different degeneracy 
parameters, they do not cross. Our spectrum s{lu) with 
u! = Arri^/2E represents the difference between the dis- 
tributions of the two flavors because the identical parts 
drop out of the oscillation equations. The only crossing 
occurs at a; = at the spectral junction between antineu- 
trinos and neutrinos. 

Flavor-dependent neutrino chemical potentials are im- 
portant inside a supernova core and also in the early 
universe if primordial neutrino asymmetries exist. The 
usual picture of synchronized oscillations and/or a syn- 
chronized MSW resonance remains unchanged. In this 
context the parametric resonance apparently does not 
lead to a novel speed-up of flavor conversions. 

Neutrinos streaming from a supernova core are not 
black-body radiation. Even if one approximates the 
fluxes of i>e , i^e and the other species by Fermi-Dirac dis- 
tributions, the effective temperatures and chemical po- 
tentials are different. Typically one finds three spectral 
crossings, the usual one at a; = and one for a nega- 
tive and one for a positive oj, see for example Ref. [32]. 
However, neutrinos streaming from a source are strongly 
anisotropic, so the parametric resonance would be sup- 
pressed by multi-angle effects as discussed in Sec. IIII Dl 
even if the spectral conditions were appropriate. 

Therefore, it is unclear if the collective motion rep- 
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resented by the parametric resonance is realized in any 
practical astrophysical or cosmological context. Perhaps 
the main insight is that even for a slowly changing or 
fixed fi the individual polarization vectors Pj need not 
follow their "single-particle Hamiltonians" even if the 
system as a whole evolves adiabatically. 

In the supernova context one is primarily concerned 
with the fate of flavor-dependent neutrino spectra as a 
function of radius, the decreasing neutrino-neutrino in- 
teraction producing a "spectral split." In the homoge- 
neous and isotropic model represented by Eq. ([T|) this oc- 
curs if /i slowly decreases so that changes from being 
/LtP-dominated to being WiB-dominated. For a single- 
crossed spectrum the occurrence of a split is explained 
by each Pi following its in the co-rotating frame. The 
spectral split occurs at a frequency Wgput corresponding 
to the final co-rotation frequency. Its value can be found 
from the conservation of B • P. In the present case of a 
double-crossed spectrum a slowly decreasing ^ leads to 
two spectral splits. Multiple-crossed spectra lead to mul- 
tiple spectral splits. The occurrence of multiple splits is 
not directly accounted for by the picture of all Pi fol- 
lowing their Hi in a single co-rotating frame and a full 
theoretical understanding of the phenomenon of multiple 
spectral splits is presently missing. 

The parametric resonance and multiple splits are not 
directly related even though this investigation was mo- 
tivated by the numerical observation of multiple splits. 
A parametric resonance for our box-spectrum requires 
the flipped middle part not to exceed half the box width 



whereas a double split occurs for any width of the flipped 
part. It is actually surprising that the large-^ oscillations 
have no apparent impact on the final sharp double-split 
that is found in the adiabatic limit of a slowly decreasing 
fi. One also finds that multiple splits are not, or at least 
not always, prevented by multi-angle effects for neutrinos 
streaming from a source. On the other hand, the collec- 
tive motion found here is "fragile" and easily suppressed 
by multi-angle effects as explained earlier. 

Perhaps the main lesson from our investigation is that 
the innocent-looking EOMs of Eq. ^ continue to sur- 
prise us with unexpected solutions, even if these solutions 
need not be of direct astrophysical relevance. Much of the 
insights gained about collective neutrino oscillations are 
owed to "numerical experiments" followed in some cases 
by theoretical explanations and even analytic solutions. 
One may well worry if our current understanding of these 
nonlinear equations is sufficiently advanced to arrive at 
robust conclusions about possible observable effects for 
example in supernova neutrino spectra. More theoretical 
work may well turn up more surprises. 
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